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Analytical Study of Twin - Jet Shielding 
Development of 3 - Dimensional Model 


I SUWARY 


Previous project reports on the analytical study of twin-jet shielding 
have dealt with the formulation of the analytical model (1) and the solution 
of the model in two dimensions (21 . The first report considered a point noise 
source impinging on a cylinder of heated flow. The solution reduces to an in- 
definite integral involving Bessel Functions. Methods to be investigated in 
order to solve the integral were outlined in the report. The second report 
discussed solution of the model consisting of a line source impinging on a 
cylinder of heated air. This two-dimensional approach eliminates the integral 
form of the solution. The trend of the estimated variation in shielding with 
frequency and azimuthal location was found to agree with the experimental data. 

In the present report, the solution for a point source impinging on a 
cylinder of heated flow is presented. The indefinite integral is solved 
approximately using a saddle of point method. Comparison of the three-dimen- 
sional model to the previously obtained two-dimensional model, indicate that the 
approximate solution of the integral is valid. 

The model is analysed in order to differentiate among the mechanims of 
shielding. The zones in which diffraction and transmission dominate are iden- 
tified. 

The model is found to compare to experimental shielding results. The 
major discrepancy between the model and experiment is felt to arise from the 
difference between the single point noise source of the model and the distributed 
source of the jet. Refinement of this noise source to a more realistic 
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representation of the real jet is one of the goals of the proposed continuation 
of this project (6). 

Further work in the present project includes investigation of the effects 
of mach number, jet temperature and source/ shielding jet separation on shielding. 


II FORMULATION OF THE MODEL 


The derivation of the model was shown in the first progress report (1). 

A stannary of that derivation is shown below. 

The mechanisms by which shielding occurs are reflection and refraction of 
sound at the boundary between the jet and the surrounding air and by diffraction 
around the jet. 

The noise source is modelled by a stationary, discrete frequency point 
source located at (r Q ,9 q , Q^) . Th< s .ieldiug jet is a cylinder of radius a, and 
is infinite in extent along the z - axis . The temperature and flow velocity are 
uniform across the cylinder cross-section. The model is illustrated in Figure 1. 
The expression for acoustic velocity potential is written for two regions; region 
I is outside the jet, region II is within the jet. 

In region l (outside the flow- incident upon the shielding jet) 

** ■ 4*tt ’ Q 0 e * 1<jt 6 ( r - r 0 ) «(«-«,) «(Z) la) 


In region I (outside the flow- reflected from the shielding Jet) 
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In region II (inside the flow) 


Where: 


** - M 2 * zz 


31 ♦ --La 

Cl ? Zt c2 *tt 


0 


(r , 8 , o,) - location of point source 
(r, 9, z) - Location of the receiver 
u> - Source frequency 
Q q - Source strength 


lb) 


2 ) 


3 



c - sound speed 

M • mach number - (jet flew speed/ci) 
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Note: Hie subscrip, refers to conditions outside the flow (ambient) , and 1 
refers to conditio, vithin the heated jet. 

The boundary conditions at the interface between the ambient air and the 
jet are: 

1) Pressure continuity 

(p) 0 ■ (p) l at r - a 
or- 

•pgC^o " " f i^t + V *z>i at r » a 3) 


2) Continuity of the vortex sheet (2) . This condition states that the 
displacement of the medium is continuous and symmetrical at the boundary; 
r ■ a. Denoting this displacement by n * (z,t) , then: 


or- 
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at r ■ a 


(n t ) 0 - (n t ♦ VOi, at r » a 
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Time is eliminated from equations 1 and 2 by assuming: 

6(r, e , z, t) - *(r, 9, z)e" iut 


The problem is reduced tc a two -dimensional formulation by the Fourier 
transform: . i r -ik.t 
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with corr e spo n ding inverse: 



Solution of the transformed equations , inclusion of the boundary conditions 
and inverse transformations yields the equation for the acoustic velocity 
potential in the far field (r»a) : 
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Ill SOLUTION OF THE MODEL 


An approximate solution of the integral in equation S is obtained using 
the Method of Stationary Phase (3) . By this method, the solution of the integral 
I z , of the fotm: 


is: 



as z ■ 

g(a 0 ) e^ChCaoHt »/*) 
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where: 

11 a Q solves h' (a) ■ o 

2) the sign in the exponential term goes as the sign of h" (t Q ) 


In OTder to solve equation S in the manner prescribed by equatioi % 
the following transformations are made. 

Let: 


1 


/ 
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then: 


1) dk, 
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k a da 
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2 ) - Mk z ■ k^ ^ — - M / 1-a 2 j 

3) K « [ k 2 - k, 2 ]** - k a 

O 0 Z 0 

1) K 1 - [ k x 2 - k z 2 ] la » k f ( — - M i 1-a 2 ’) 2 - 1 * a 2 l ** 
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the integral in equation S : 


becomes: 


C 

1: 
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In the acoustic far field, r»l. Assuming that k 2 «k Q , argument of the 
Hankel function becomes large: 

k ira» 1 
o 

and the Hankel function is approximated by: 
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fran figure 1, it is seen that: 

r * R cos 8 and z * R sin 8 

where: 

R ■ distance from origin of the coordinate system to the receiver 
8 ■ angle between vector R and the projection of the vector R on the 
x - y plane. 


The integral is then: 


f, 


-i — (m 


irR COS 8 




R( V i-a 2 sin 8 '♦a cos 8), 



F (a) da 
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comparing equation 10 to equation 6 : 
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thus: 


then: 


and: 


z ■ k R 

o 

h(«) ■Vi-a 2 sin B ♦ a COS 6 
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a is evaluated from: 
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h' (a Q ) ■ cos S 
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and: 






cos B ' 
sin 2 B 


F m (c° s B ) 


The integral, as k Q R >> 1, becomes: 
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From which the expression for the far field acoustic velocity potential 
is evaluated: 

Qq© iut • -imir ik R 
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Where: 
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invert: 


Ti ■ j- * M sin 8 
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The total acoustic velocity potential in the far field consists of two parts: 

1. The velocity potential of the incident wave: 

0e' iut . -inn ik 0 K 
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2. 'the velocity potential of the scattered wave: 
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| terms in equation 14) within brackets j 

As a check; when r Q • o (source is located at the origin) , the velocity 
potential of the incident wave reduces to: 


16 ) 



sines; 
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This is the expression for the velocity potential in the far field for 
radiation from a point source of strength, located at the origin of the co- 
ordinate system. 


IV RESULTS 


The total sound pressure is evaluated from the acoustic velocity potential 


by: 


P T' -»„& T 


The incident sound pressure is evaluated from equation IS: 
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The normalized sound pressure is formed from the magnitude of the ratio of 
the total sound pressure to the iixident sound pressure. This normalized 
sound pressure is a measure of the influence of the shielding jet on the noise 
source directionality. A value of normalized sound pressure less that ore in- 
dicates noise reduction; and a value of normalized sound pressure greater than 
one indicated amplifications. 

1. Comparison of 3 - dimensional Model co 2 - dimensional Mod;' 


The 3 - dimensional model is compered to the 2 - dimensional model developed 
in a previous report (2). The conditions imposed are: 

a) M ■ Mach nunber ■ 0.0 

b) B ■ 0.0 

C) 0-0 ■ w 

' 0 

Under these conditions, the receiver is located in the x-y plane of the 
source, on the side of the jet opposite the source. It is expected that the 
normalized sound pressure reduces to the two dimensional case investigated in 
a previous progress report (2). The results obtained are shown in figure 2 
and are compared to plots obtained previously. The plots of normalized sound 
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pressure against non-dimensional frequency parameter , k Q a, are comparable from 
the two studies. This indicates that the approximation of the integral employed 
in the present study is valid. 

2. Variation of Normalized Sound Pressure with Angle to Jet Axis. 

In order to differentiate among the mechanisms of shielding , the normalized 
sound pressure is evaluated at low frequencies under the following operating 
conditions. 

a) Vj * Jet velocity = 1552 feet/sec 

b) Tj = Jet temperature = 1238°R 

c) T q = Ambient temperature * 530° R 

d) S/D = Spacing parameter = ratio of distance from source to center of 

jet and jet diameter. =2.67 

The normalized sound pressure is evaluated in the shielded zone of the jet. 
Thus, 9 -Q q = tt and o* 8 < ir/ 2 . For the purposes of comparison to other data, 
the coordinates are shifted to a set centered in the source, and illustrated 
in figure 3. In the nomenclature adopted in the figure, 4» n - o° is on the 
z-axis of the noise source, and $ = 90° is on the side of the jet directly opposite 
the source. 

The contours of normalized pressure for 10° ^ ^ < 90° are plotted against 

non-dimensional frequency parameter in figure 4. The non-dimensional frequency 

k a is varied from 0.1 to 1.0. 
o 

At angles close to the z-axis of the noise source, 10° < < 40°, the effect 

of increasing the angle is to shift the normalized sound pressure curve to 
lower frequencies. As the angle increases, the line of sight distance from the 
source to the jet becomes smaller. From barrier theory arguments developed in 
a previous report (2), the expected result of decreasing the spacing between the 
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source and the shielding jet is to increase the attenuation of the lower frequency 
sounds. Thus, in the range 10 9 < < 40°, diffraction around the jet dominates . 

In the range of incident angles 40° < # n < 90°, the trend reverses, and in- 
creasing the angle of incidence shifts the curves toward higher frequencies. 

In this range, the sound transmitted through the jet is combining with the dif- 
fracted sound. This conclusion is consistent with analysis by Yeh (5). Yeh 
considered a sound ray incident upon a uniform sheet of thickness, D, in which 
the flow and temperature are constant. The sound ray is inclined at an angle 
to the axis of the flow sheet. The model is illustrated in figure 5 . The 
temperature difference and flow cause the incident sound ray to be refracted, with 
the angle in the jet, <|>j, evaluated from: 

cos - C J cos *n 

c o - V J cos *n 


Transmission through the jet is minimized when the refracted ray becomes 
parallel to the jet or <j»j * o. The critical value of the incident ray, i|/ nc , 
corresponding to t|/ T * o is: 

C o 

COS lb 


nc 


r + v 

J J V J 


The jet parameters: 

C Q - 1128 feet/sec 
Cj * 1725 feet/sec 
Vj ■ 1552 feet/sec 

result in: 
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Yeh's model estimates that all the incident sound in the range i 

T n nc 

is transmitted through the jet, while no transmission occurs for sound incident 
at an angle * n < * nc . 

Referring to figure 4, the normalized sound pressure decreases more rapidly 
in the range of angles 40° < f < 70° than it does in the range 70° < < 90°. 

This indicates that sound incident at i|» n > 70° is transmitted through the jet; 
and that the contributions from the transmitted sound decreases rapidly for 
^n < 70 °' At incidence, < 40°, the diffracted sound contribution 

dominates . 

3. Comparison to Experimental Results. 

The analytical model is compared to experimental results by Kantola (4) 

and plotted in figure 6. The condition evaluated is: 

a. Tj = 1238° R 

b. Vj * 1519 feet/sec 

c. S/D = 2.67 

d. * n - 30°, 60° 

The analytical model shows some agreement in foim with the experimental data. 
The onset o. hielding occurs at lower frequency for * 30° than it does for 

<|> n * 60°. As the jet axis is approached, the agreement between the analytical 
and experimental improves, with the significant difference being at high fre- 
quencies. The analytical model shows that the normalized sound pressure increases 

for k Q a > 10. This indicates a upper frequency limit on shielding. The 
experimental results indicate no such upper limit in the range of frequencies 
investigated. 

The major discrepancy between the analytical and experimental results occurs 
at greater angles of incidence. This is felt to be due to the point source 
approximation of the source jet. The jet noise source is more realistically 
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considered to be a distribution of point sources along the jet axis. The 
downstream sources are contributing to the sound pressure measured in the 
shielded zone. From the previous discussion on transmission though the 
shielding jet, the sound incident on the jet at ij» n > 70° is transmitted more 
readily than sound incident at ^ < 70°. In general then, sound is transmitted 
upstream more efficiently than downstream because the angle of incidence is 
greater. Thus, the receiver at ip n ■ 60° "sees" more of the downstream noise 
sources than does the receiver located closer to the jet axis. The point source 
approximation for the source jet breaks down as the angle of incidence increases. 
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V. CONCLUSIONS 


The 3 - dimensional model of twin jet noise shielding has been obtained for 
a point noise source impinging on a cylinder of heated flow. The validity of 
mathematical approximation required in the solution is tested by comparison of 
the 3 - dimensional model to a previously obtained 2 - dimensional fbnnulation. 

In the shielded zone, the normalized sound pressure curve is shifted toward 
higher frequencies as the angle of incidence from the jet axis is increased. 

This trend is expected from barrier theory. In this range of angles, diffraction 
around the jet dominates. For the jet parameters chosen, diffraction is the 
dominant mechanism for <i» n < 40°. As the angle of incidence increases further, 
transmission of sound through the jet contributes to the diffracted sound. At 
angles of incidence greater than a critical value defined by: 


cos f 


nc 


c j + v j 


Transmission through the jet is the dominant mechanism. 

The analytical model compares to experimental results more favorably for 
angles of incidence close to the jet axis. The discrepencies that occur at 
greater angles of incidence are felt to arise from the difference between the 
single point source approximation in the model and the distributed source con- 
figuration of the real jet. 


- 16 - 



REFERENCES 


L. Gerhold, C. ; 'Development of an Analytical Model of Twin Jet Shielding". 
Progress report or NASA Grant NAG-1 -11. June 18, 1980 

2. Gerhold, C. ; "Analytical Study of Twin Jet Shielding - Two-Dimensional 
Model". Progress report for NASA Grant NAG-1-11. August 14, 1980 

3. Jefferies, H. and Jefferies , B. ; Methods of Mathematical Physics . 
Cambridge University Press, 1956 

4. Kantola, R.A. "Shielding Aspects of Heated Twin Jet Noise", AIAA 4th 
Aeroacoustics Conference, Atlanta, Georgis, October 3-5, 1977 

5. Yeh,C.; "A Further Note on the Reflection and Transmission of Sound 
Waves by a Moving Fluid Layer". Journal Acoustical Society of America 
vol. 43, no. 6, pp. 1454-1455, 1960 

6. Gerhold ,C. ; "Analytical Study of Twin -Jet Shielding Progress to Date and 
Proposed Project Continuation". Progress report for NASA Grant NAG-1-11. 
October 1,1980 


L 


Receiver 

(r,0,z) 





- 18 




Figure 1 Source - Receiver Geometry for 
Jet Shielding Analysis, 
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igurc 3 Definition of Angle in Coordinates centered on Noise souro^ 
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Figure 5. Sound Transmission Through Moving Heated Layer 
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